Model independent tests of the standard cosmological model 



Arman Shafieloo 1 & Chris Clarkson 2 

1 Department of Physics, University of Oxford, 1 Keble Road, Oxford, OX13NP, UK. 

2 Centre for Astrophysics, Cosmology and Gravitation, and, Department Mathematics 
and Applied Mathematics, University of Cape Town, Rondebosch 7701, South Africa. 

The dark energy problem has led to speculation that not only may ACDM be wrong, but that 
the FLRW models themselves may not even provide the correct family of background models. We 
discuss how direct measurements of H(z) can be used to formulate tests of the standard paradigm 
in cosmology. On their own, such measurements can be used to test for deviations from flat ACDM. 
When combined with supernovae distances, Hubble rate measurements provide a test of the Coper- 
nican principle and the homogeneity assumption of the standard model, which is independent of 
dark energy or metric based theory of gravity. A modification of this test also provides a model 
independent observable for flatness which decorrelates curvature determination from dark energy. 
We investigate these tests using Hubble rate measurements from age data, as well as from a Hubble 
rate inferred from recent measurements of the baryon acoustic oscillations. While the current data 
is too weak to say anything significant, these tests are exciting prospects for the future. 



1. INTRODUCTION 

There is a growing realisation in cosmology that the 
dark energy problem means that we have to test the 
fundamentals of our cosmological model as rigorously as 
we can. Despite the fact that the standard flat ACDM 
paradigm works so well on many fronts, the unknowns 
at the heart of it have led us into the territory of trying 
to invert observations to ascertain directly any temporal 
evolution of dark energy - whether it is in fact modi- 
fied gravity, global inhomogeneity such as a void model, 
quintessence and so on. Normally in physics, we would 
hope to work the other way round: a model based on 
known physics can predict outcomes of an experiment 
which can then potentially falsify the model (or theory). 
Unfortunately, once we move away from the comfort of 
ACDM we have little physics to guide us at all; worse, 
many suggestions appear degenerate or even unfalsifiable. 
Generalised gravity theories, quintessence models, void 
models, etc. all have free functional degrees of freedom, 
and not just free parameters. So dark energy reconstruc- 
tion tries to deduce functional degrees of freedom directly 
from data sets we have. 

In particular, if Hubblc-scale inhomogeneity is respon- 
sible for the apparent acceleration, then all bets are off as 
to what the effective equation of state 1 could be. In the 
simplest void models which describe this [TJ[2] , a violation 
of the Copernican principle - that we are not located at 
a special place in the universe - is required, because they 
place us very near the centre of symmetry. Clearly such 
a radical suggestion warrants much more study; for ex- 
ample, it is not yet understood how perturbations evolve, 
although the generalised Bardeen equation is known [3]. 



One important aspect of the reconstruction approach 
will rely on formulating consistency tests of the standard 
paradigm. Ideally these tests would be constructed so 
that signals of deviations from (flat) ACDM are easy to 
spot and quantify. If it's possible, we would like these 
tests to be independent of the standard paradigm, and 
rely as little as possible on assumptions that we normally 
take for granted. For example, consistency between the 
background evolution and the evolution of perturbations 
can be used to rule out whole classes of quintessence mod- 
els [3]. Using perturbations in this way is only possible 
within given classes of dark energy, however, because of 
assumptions about the sound speed which must be made, 
amongst other things. Tests which may be formulated 
purely from background observables would not suffer this 
difficulty. 

Here our principle aim is to investigate the curvature 
test, introduced in [5]. This is a test which uses back- 
ground observables only, and tests for consistency of the 
FLRW paradigm itself. In FLRW models, the luminosity 
distance may be written as 



d L (z) = C(1 ^l_ sin ( J^h [" dz' 
M ; iW-^fc V Jo H(z') 



(1.1) 



where Qk is the curvature parameter today, and the ex- 
pansion rate H(z) takes the value H = H(0) today. 
(Note that this is actually valid for either sign of curva- 
ture; for ftk = we can take the limit fi^ — > 0.) The area 
distance is defined using cLl, = (1 + z) 2 dA, and another 
distance measure we will use is D = (1 + z)dA/{c/ Hq), 
which is the dimensionless comoving distance. We may 
rearrange Eq. (1.1) to give an expression for the curva- 



ture parameter in terms of H{z) and D(z) [5]: 



1 That is, the equation of state which reproduces some particular 
family of observations, such as the distance modulus, if we utilise 
a FLRW model. 



[h(z)D'(z)Y 



(1.2) 



where h(z) = H(z)/H and ' = /z. This tells us how to 
measure the value of the curvature parameter today from 
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distance and Hubble rate observations, independently of 
any other model parameters or dark energy model or field 
equations. Remarkably this tells us the curvature today 
from these measurements at any single redshift, provided 
that D(z) and H{z) come from unrelated families of ob- 
servations. If the underlying geometry of the background 
cosmology is FLRW, &k(z) will be measured as constant 
and so be independent of the redshift of measurement. 
As a consistency test of the standard paradigm it has all 
the right ingredients. It is a straightforward comparison 
of different obscrvablcs, both of which can be formulated 
fairly independently of contamination from the standard 
model. At the simplest level, then, we can use this to 
see if our distance measurements are consistent with our 
H(z) measurements, if we choose to stick within the stan- 
dard paradigm. But we can get more information than 
that. 

For measurements of &k{z) not to be constant as 
a function of redshift we are looking at pretty radical 
changes to the standard model - an alteration to the un- 
derlying homogeneous and isotropic FLRW models them- 
selves. One important situation where this will happen 
is in void or Hubble-bubble models. These models fit 
the SNIa distance modulus by having us near the centre 
of a large Hubble-scale depression in the energy density 
and accompanying curvature. These models violate the 
Copernican principle which lies at the heart of cosmol- 
ogy. If isotropy about us were determined at all red- 
shifts, this curvature test may then be used as a direct 
test of the Copernican principle. However, deviations 
from &k{z) =const. will be signalled in any non-FLRW 
cosmology; in particular, in Swiss-cheese models which 
satisfy the Copernican principle on very large scales 
will have an &k{z) which oscillates with redshift. Fur- 
thermore, backreaction and averaging will induce non- 
constant 0k{z) [6]. (Note that in all these models it 
won't necessarily be saying anything specific about the 
curvature.) 

A similar background test which is specifically to test 
any deviations the flat ACDM scenario was introduced 
independently in [7J|H]- The Om diagnostic defined by [7J: 



Om(z) = 



h 2 (z) - 1 _ 

(1 + z)3 - 1 ~ [(1 + z)3 - l]D'{zf 



1-D'(z) 2 



(1.3) 



is a constant at different redshifts in a spatially flat 
ACDM model, equal to today's value of the matter den- 
sity parameter f2 m . In a similar spirit to f2j- we can mea- 
sure the variables on the rhs, and expect our measure- 
ments to yield the same number regardless of redshift. 
Fig. [T] shows the kind of behaviour Om(z) exhibits from 
different FLRW models. From this we can also derive [5] : 



(z) = 2[(l + zf-l]D"(z) 

+3(l + z) 2 D'(z)[l-D'(z) 2 ] 
= for all flat ACDM models. 



(1.4) 




redshift z 

FIG. 1: A plot of Om(z) in different models, obtained us- 
ing distances. The figure shows the range of behaviour from 
curvature in each fan of curves (key, top left). The three 
grey fans show how differing £l m values interact with curva- 
ture for ACDM (key, top right). The effect of changing w by 
a constant is illustrated in the red and brown fans. Similar 
curves are found if we use h rather than D to form Om(z), 
but without the cross-over at high redshift. It is clear that, 
although Om(z) deviates from a constant because of curva- 
ture this deviation is not strong for small curvature values; 
therefore, Om(z)'s utility lies in picking up non-A behaviour. 



The main issue in using the curvature test lies in ob- 
serving H(z) independently of distance measurements. 
For example, H{z) is often reconstructed from super- 
novae data [5], but this implicitly uses Eq. (1.2 1 and so 
can't be used here. There are other methods, two of 
which we use here. Measurements of the Baryon Acous- 
tic Oscillations give the 'volume distance' which is 



D v (zf = 



c 

Ho 



zD(z) 2 



(1.5) 



Testing for deviations from Jz?(z) 
nate test of flat ACDM. 



provides an alter- 



from which we can derive H(z) if we have an another 
method to determine D(z). Alternatively, we can refor- 
mulate the curvature test directly in terms of Dy{z) and 
D(z). This method is not entirely independent of cosmo- 
logical model, however, because it assumes that there is 
no scale dependence in the evolution of perturbations - 
this is not the case in inhomogeneous models where the 
curvature is varying spatially [3]. 

Spectroscopic age dating of passively evolving galaxies 
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give another measure of H(z) [3T], which we use here. 
This method measures 



1 



(l + z)H(z) 



Other possible methods include measurements of the 
time drift of redshifts [TU] and the dipolc in the SN ob- 
servations [IT], but these are not considered further here. 
For distance measurements we shall use the latest com- 
plete SNIa data set [15] . By using the smoothing method 
presented in [T2TH4] we can turn this data into D(z) and 
D'(z) without relying on any hducial cosmological model 
or parameterisation. We use &k{z) = fife in Eq. (1.2) 



and Om(z) in Eq. (1.3) as two diagnostics which we at- 



tempt to reconstruct using these independent datasets. 
Any deviations from constant will signify problems with 
the FLRW models themselves in the case of &k{z), and 
problems with flat ACDM in the case of Om(z). Finally 
we make use of a third test: ^(z)Z?(z) 2 may be used 
as a model independent test for flatness if we look for 
deviations from zero. 



2. CONSTRUCTION OF THE DIAGNOSTICS 
A. Finding D(z) and D'(z) from supernovae 

To find D(z) and D'(z), we use the smoothing method 
implemented and applied to supernovae data by |12[ 1 1 3] . 
In addition we have done some modifications to the 
smoothing method to make it error-sensitive that results 
to a better fit to the data [TT]. This method is a com- 
pletely model independent approach to derive the di,(z) 
relation directly from the data, without any assumptions 
other than the introduction of a smoothing scale. The 
only parameter used in the smoothing method is the 
smoothing width A, which is constrained only by the 
quality and quantity of the data, and has nothing to do 
with any cosmological model. The smoothing method is 
an iterative procedure with each iteration typically giv- 
ing a better fit to the data. It has been shown in [T2 [T3] 
that the final reconstructed results are independent of 
the assumed initial guess, d^Zi) 3 below. 

The modified smoothing method (error-sensitive) can 
be summarized by the following equation [14 : 



lnd L {z,Af = ln d L {z) 9 

+n(z) £ d _iMn exp 



hr 



2 / 1+Zj 
1+z 



2A 2 



where 
N{zY 1 = ^]exp 



lrr 



2A 2 



The absolute brightness of the supernovae is degener- 
ate with Hq since the observed quantity is the distance 
modulus fJ-(z). The outcome of the smoothing method 



(1-6) is therefore H d L {z)/c = rf L GC (z) = (1 + z)D(z 



We 

should note that the actual value of Hq is quite impor- 
tant in our analysis; supernovae surveys usually assume 
the value of Hq in their data which can bias our results. 
In this paper we use the recent 'constitution' supernovae 
data [H] that include 397 data points and we choose 
A = 0.30. This value of A is half of the value used in 
[T3] in analysis of the SNLS data [IS] , simply because we 
have almost four times as many data points in the consti- 
tution sample in comparison with SNLS data. Complete 
explanation of the relations between the A, the number 
of data points, quality of the data and the reconstructed 
results can be found in [T2J H3] • Our best reconstructed 
result using smoothing method gives \ 2 = 459.4 to the 
constitution data, which is a significantly better fit than 
the best fit x 2 derived by using the CPL parameteri- 
zation w(z) = w + wijfi [IS], with Ax 2 « 1.7 [17] , 
By using our smoothing method we reconstruct a set of 
d T £°(z) in the redshift range of the supernovae data. We 
consider only the reconstructed results that have a bet- 
ter fit than the best fit spatially flat ACDM model with 
X 2 = 465.63 (fiom = 0.288) [17]. We show these curves 
in Fig. [2j We should note that though all curves we 
have shown in Fig. [2] have the \ 2 better than the best fit 
ACDM model, but we should consider that the A% 2 be- 
tween these reconstructions of the expansion history and 
the best fit ACDM model is not that large to conclude 
that flat ACDM has a strong discrepancy with the data 
(maximum A\ 2 is around 6). 

Once we have reconstructed D(z) finding D'(z) is a 
matter of differentiating our D(z) curves. Because the 
procedure is iterative the errors on D'{z) may be esti- 
mated from the set of D(z) curves which we consider 
with x 2 < 465.63. 



B. Finding h(z) from the BAO 

1. D v {z) from the BAO and CMB 

From direct observations, there are two CMB param- 
eters, R, the scaled distance to recombination, and £ a , 
the angular scale of the sound horizon at recombination, 
with measured values that are nearly uncorrelated with 
each other [19]. These parameters are defined as: 



where g?l(z) is the data, N(z) is the normalization factor, 
dh^i) 9 is the initial guess model and A is the width of 
smoothing. 



and 



(2.1) 



R= \l VL m Hlr{z C MB) 

^r(z C MB) 
r s {zcMB) ' 



(2.2) 



(2.3) 



where r(z) = cD(z)/H is the comoving distance, and 
zcmb = 1090. It is also possible to put constraints on 



4 



x 



o 

a 

a 
c 



o.io- 



0.05- 



-0.05 



-0.10 




2. Derivation of h(z) 



redshift z 

FIG. 2: The best fit curves to the constitution supernovae, 
shown here with vertical error bars reduced by a factor of 10 
for clarity. The distance modulus is shown with the best fit 
flat ACDM subtracted off (with H = 65 kms^Mpc" 1 ). The 
curves all have a better fit to the data than ACDM, shown 
by the dashed horizontal line. Note the lack of SNIa above 
z ~ 1 means that D(z) is not well constrained there. 



t s{ z cmb)) the comoving sound horizon at recombination, 
directly from CMB data independently of the assumption 
of the cosmic curvature or model of dark energy [TU] : 



r s (z C MB) 



c 



dz 



h(z)> 



(2.4) 



From D 



and d r £° we may estimate h{z) from 



h(z) 



c 



zd r £ c 



(zf 



(l + zfD v {zY 



(2.6) 



To calculate h(z), therefore, we also need to know the 
value of Hq. The uncertainties of Hq can affect the re- 
construction of the cosmological quantities. As we see in 
eqs. 1.2 and 1.3, to test cosmological models we need to 
know h(z) = H(z)/Hq. So only knowledge of H(z) is not 
enough and the value and error-bars of H are also very 
important and they affect the reconstructed h(z) and its 
error-bars. In fact the first term in the right hand side of 
eq. 2.7, is the error from the uncertainties of Ho- For our 
analysis we choose different constraints of H Q from dif- 
ferent observations. For each reconstructed set of d r £ c (z) 
from supernovae data we can derive o'htz) '■ 



&h(z) 




a H 



dh(z) 
dD v (z) 



a D v (z) 
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D v (z)_ 





From <Jh( z )i we can evaluate <jg k and <Jom'- 



and 



2d r £ c '(z) 2 h{z) 

2h(z) 
(l + z) 3 -l 



'h(z) 



>h(z)- 



(2.7) 



(2.8) 



(2.9) 



Note that these errors are for a given d^ c (z); errors re- 
sulting from the SN la measurements are found by cal- 
culating h(z) for many d r £ c (z). 



where c s is the sound speed and the derived value of 
r s{ z cMB) from WMAP data is r s (z C MB) = 148.55 ± 
2.60Mpc Ql]. 

The BAO give a measurement of &{z) = rs{ ^^f ) at 
two redshifts of z = 0.35 and z = 0.20 gD]: ^^o) = 
0.1905 ± 0.0061 and r -f^W^ = 0-1097 ± 0.0036. Com- 

Dy (0.35) 

bining these measurements with the value of t s (zcmb) = 
148.55 ± 2. 60Mpc from WMAP data we can estimate the 
values of Dy(0.20) and £>y(0.35). Since we are dealing 
with two completely different observations we can assume 
that they are almost uncorrelated. Hence we can derive 

a D v (z)- 



&D v {z 



dDyjz) 

dr s (z CMB ) 



a r 3 {z C MB) ' 



dD v (z) 



&m(z) 



d@{z) 

(2.5) 

which results in L> y (0.20) = 779.79 ± 28.45Mpc and 
LV(0.35) = 1354.14 ± 50.36Mpc. 



C. Finding h(z) from age data 

In gTJ [52] H{z) is derived by using the relative ages 
of passively evolving galaxies. We can simply derive h(z) 
and c^( z ) using this data: 



(2.10) 



This assumes that our measurements of Ho and H(z) are 
uncorrelated. We can now estimate Om{z) and Gk{z) in 
a same way as we did in the previous section by using 
D(z) and D'(z) from supernovae data. 



3. RECONSTRUCTED DIAGNOSTICS 

A. Om(z) and J£(z) from supernovae distances 

The tests for flat ACDM discussed above may be 
utilised directly using distance data and its derivatives, 
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redshift z 

FIG. 3: Tests for flat ACDM, Om(z) and Jzf(z). The smoothing method allows us to construct these directly. The curves 
shown are all a better fit to the data than flat ACDM, which is shown as the dashed line in the plots. These two plots are 
made from the same reconstruction. We have shown Jz?(z) up to z = 1 because the error-bars at the higher redshifts are so 
large. Jz?(z) is derived by using the second derivative of the results from the smoothing method were Om(z) is using the first 
derivative, so the error-bars for Ji?(z) are larger. 



which we show in Fig. [3] Both Om{z) and J£{z) tell 
us useful information about possible deviations from flat 
ACDM. For example, at z « 0.3, Om(z) appears con- 
sistent with f2 TO sa 0.3; ££{z) tells a different story as it 
deviates from zero there. For both curves there is a pref- 
erential deviation from the concordance model. While 
increasing behaviour of Om(z) at low redshifts suggests 
larger w(z) for dark energy in this range, the large er- 
rors at high red shifts above z ~ 1 indicates poor quality 
of the data. While these curves are preferred over flat 
ACDM, they are not strongly preferred enough to get too 
excited about. It has been discussed before that even by 
using parametric methods in the reconstruction of the 
expansion history, the Om diagnostic would not be very 
sensitive to the choice of parameterization; as we see here 
our results for Om(z) using a non-parametric smoothing 
method are very similar to reconstructed Om(z) using 
a series ansatz for w(z) [T7j. It is also rather surpris- 
ing that the curves derived here, which are model inde- 
pendent constructions of Om{z), have similar features to 
the corresponding curves we find in void models with a 
smooth centre 0. 



B. The curvature test and Om(z) using Hubble rate 
measurements 

The diagnostics are sensitive to the value of H . 
Assuming free priors on curvature and also equa- 
tion of state of dark energy, Ref. [23] find H = 
68.7 ± 2.0 kms _1 Mpc _1 from a combination of 
WMAP5+BAO+high-z SNe data. The most recent con- 
straints on Hq from the observation of 240 long period 



Cepheids by Ref. [35] give us instead H = 74.2 ± 3.6 
kms _1 Mpc _1 . The latter of these results is more suit- 
able for us because it makes fewer assumptions on the 
underlying cosmology; we will consider @k(z) and Om(z) 
using both of these measurements however. 

In Fig.|4]we present the Om(z) diagnostic for h(z) data 
from the BAO and age data (note that age data doesn't 
use the SN distances, but that the BAO does, when we 
reconstruct h(z)). These data points are essentially a 
remapping of the data through the function Om{z). Note 
that the errors shrink with distance, and that the data 
are clearly consistent with Ora(z) =const. 

In Fig. [5] we show the reconstructed {?k{z) using both 
BAO and age data to evaluate h(z) using both values of 
Ho (centre). For low redshift the errors are very large 
indeed, but for high redshift this test shows considerable 
promise. In the zoomed picture on the right, we can 
see that &k(z), while not tightly constrained by current 
data, the error bars are reasonable in comparison with 
Om(z), as most of the extra uncertainty comes from er- 
rors in the reconstructed D'(z). We can see by eye that 
0k[z) =const. is consistent with the data, but that neg- 
ative values are slightly favoured. This is actually re- 
flection of the fact that Om(z) deviates strongly from 
constant at high redshift, which in turn reflects the fact 
that the SNIa data are not sufficient there. 

The errors at low redshift indicate that &k{z) =const. 
will not be a useful test in this regime. However, the 
numerator of ffk(z), gives an important measure of devi- 
ations from flatness, irrespective of the dark energy equa- 
tion of state or theory of gravity (see also [5]). In the left 
panel of Fig. [5] we see that for moderate redshift data 
points we have the seed of a model independent test for 
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FIG. 4: Resultant Om(z) using two values of Ho- The points 
indicated with an arrow show the SN (constitution) data com- 
bined with the recent BAO (SDSS-DR7) and CMB (WMAP) 
data. The other points show the same using age data. The 
thick vertical lines represent the spread in the constructed 
point arising from uncertainty in the reconstructed distances; 
the thin vertical lines represent the errors arising from H(z) 
and Ho- 



flatness which will become a robust test as the data im- 
proves. In general, measuring sign(^(z)D(z) 2 ) will give 
a robust and model independent determination of the 
sign of the spatial curvature. 

The reason for the inefficiency of the Copernican test 
at low redshift may be understood as follows. Given two 
families of independent observations, {D(z), H(z)}, say, 
we may reconstruct separate FLRW models with their 
own parameters from each family of observations, modulo 
all the usual degeneracies. We may think of &k{z) as the 
mismatch between these reconstructed models. For small 
redshifts, this implies that, 



(D) 



9. 



(H) 



l {D) w {D) 

DE 



n (H) w {H) 



-O(z 



(3.1) 



where ftnF^ is the curvature parameter today evaluated 
from distance measurements, etc. The next term in the 
series is a constant and gives the value of the curvature 
parameter today when the two observables give the same 
model. For small redshift then, any errors are diver- 
gent. This limits the utility of this particular diagnos- 
tic to moderate to large redshift, where we have seen it 
shows potential. More generally, it will be useful to use 
the fact that 



<#(z) = l + h 2 (DD" - D' 



hh'DD' 



(3.2) 



must be zero in any FLRW model [5]. This doesn't di- 
verge at low redshifts; rather we have 



(D) 
% 



% 



z + 0(z 2 



(3.3) 



which is much more stable. This comes at the expense 
of requiring derivatives of Hubble data, which we can 
imagine being dealt with using the smoothing method 
we have used here for D(z) once the data becomes good 
enough. 



4. DISCUSSION AND CONCLUSIONS 

We have investigated some model independent con- 
sistency tests of the FLRW models using a mixture of 
the latest complete SNIa data and some Hubble rate 
data. The SNIa data is turned into the dimensionless dis- 
tance function D(z) using a model independent smooth- 
ing method. One source of the H(z) data is from the 
BAO measurements of Dv{z), which relies on assump- 
tions about how perturbations evolve in the standard 
FLRW paradigm. The other is from the relative ages of 
passively evolving galaxies, which has far fewer assump- 
tions about the underlying cosmological model built in 
to it. 

We have considered three tests: 

• Tests for flat ACDM: Om(z) and 3?(z). The first 
of these was investigated using both Hubble rate 
data and distance data independently. When us- 
ing Hubble rate data, everything is consistent with 
flat ACDM. Using distance data, there seems to be 
slight evidence for deviation from flat ACDM. We 
should note that a larger Om(z) at low red shifts, 
as we see in Fig. |4j is indicative of a larger w(z) 
for dark energy (considering the effect to be from 
dark energy, and not due to curvature) and suggest 
that some decaying models of dark energy (to dark 
matter or something else) must have a good fit to 
the data. This is in agreement with [17] . However 
larger reconstructed Om(z) at high redshifts (de- 
spite of the fact that the quality of the data is poor 
in this range) seems to be related to the apparent 
extra brightness of the supernovae data (in compar- 
ison with ACDM model) at z > 1 as reported and 
studied in [23] , The second test Jzf (z) relies on sec- 
ond derivatives of distance data, but is surprisingly 
useful even with present data; here, deviations from 
flat ACDM are seen with mild significance. It is in- 
teresting to note the close resemblance of the curves 
we have reconstructed here in a model independent 
way, to those same functions in void models [5] . 

• Test for the sign of curvature: &k{z)D{z) 2 is a test 
for flatness, which is independent of dark energy. 
It requires two independent data sets, and seems 
most useful at the moment for z < 1. For z > I 
the errors on SNIa distances dominate. There is no 
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FIG. 5: Resultant &k(z) using two values of Ho. The points indicated with an arrow show the SN (constitution) data combined 
with the recent BAO (SDSS-DR7) and CMB (WMAP) data. The other points show the same using age data. The thick vertical 
lines represent the spread in the constructed point arising from uncertainty in the reconstructed distances; the thin vertical 
lines represent the errors arising from H(z) and Ho- The faded points are two age data points which are at such high redshift 
they have no corresponding SNIa, so should not be taken too seriously. 



evidence for deviation from flatness, though this is 
not surprising because it relies on Hubble rate data. 

• Copernican test for homogeneity: &k(z) is a model 
independent test for the FLRW models themselves, 
which is independent of dark energy or theory of 
gravity. It requires two independent data sets, and 
is very promising for z > 1. Again, because it re- 
lies on Hubble rate data, there is no evidence for 
deviations from FLRW, nor indeed from flatness. 
Note that for 0.8 < z < 1.6 the data shows a mild 
trend towards closed models. This is likely of no 
significance. 

It is early days for tests as sensitive as these because 
they require non-parametric construction, and we are still 
dealing with only a few hundred data points. We can eas- 
ily envisage reaching thousands of data points in the next 
few years at which point these tests will come into their 
own as foundational tests for the concordance model. If 
the concordance model is correct then these tests will 



provide a useful way to quantify how confident we are in 
a more absolute way than we can at present. If, on the 
other hand, we find deviations for any of these tests, then 
they will help point us in other directions. In particular, 
if the Copernican tests Gk(z) and ^{z) fail then we have 
a serious problem for cosmology, and the standard FLRW 
models. But with that would come a strong indication 
that dark energy cannot be due to modified field equa- 
tions on either side; rather, we would know that there is 
something wrong with the spacetime metric we are using 
or that we are not describing the universe in the mathe- 
matically correct way. 
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